Gap maps and intrinsic diffraction losses in one-dimensional photonic crystal slabs 
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A theoretical study of photonic bands for one- dimensional (ID) lattices embedded in planar 
waveguides with strong refractive index contrast is presented. The approach relies on expanding 
the electromagnetic field on the basis of guided modes of an effective waveguide, and on treating 
the coupling to radiative modes by perturbation theory. Photonic mode dispersion, gap maps, 
and intrinsic diffraction losses of quasi-guided modes are calculated for the case of self-standing 
membranes as well as for Silicon-on-Insulator structures. Photonic band gaps in a waveguide are 
found to depend strongly on the core thickness and on polarization, so that the gaps for transverse 
electric and transverse magnetic modes most often do not overlap. Radiative losses of quasi-guided 
modes above the light line depend in a nontrivial way on structure parameters, mode index and 
wavevector. The results of this study may be useful for the design of integrated ID photonic 
structures with low radiative losses. 

PACS numbers: 42.25.Fx, 42.70.Qs, 42.79.Dj, 42.82.Et 



I. INTRODUCTION 

Photonic crystals embedded in planar dielectric waveg- 
uides, also known as photonic crystal slabs, are inten- 
sively investigated as a promising route for the tailoring 
of photonic states. 1-40 Indeed, propagation of light can 
be controlled in these systems by the dielectric disconti- 
nuity of the slab waveguide in the vertical (z) direction 
and by the photonic pattern in the xy plane. The geom- 
etry of a patterned waveguide gives considerable freedom 
in designing photonic structures (periodic or containing 
defects) that can be realized at near-infrared or optical 
wavelength by lithography and etching. 

Most experimental investigations of photonic crys- 
tal slabs with a two-dimensional (2D) or one- 
dimensional (ID) periodic lattice concern in-plane 
transmission^^S^ or surface reflectance/transmittance 
measurements^i^SSi^SiMiS with the purpose of deter- 
mining the photonic gaps and the band dispersion. Struc- 
tures containing defect states like linear waveguides in 
2D lattices^ 6 - or microcavities in ID systems^ 2 * 3 ^ . are 
also being investigated. On the theoretical side, the 
study of photonic crystal slabs has been undertaken 



with plane-wave expansion, 
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methods, 9-25 29 finite-difference time-domain (FDTD) 
calculations ) 15 ' 22 ^?' 28 modal methods ; 12 ' 1 ?' 18 ' 19 and per- 
turbative approaches Recently, a finite-basis ex- 
pansion method has been introduced. 36-38 Most of these 
papers concern 2D structures, either periodic or with 
linear defects. The theoretical study of ID structures 
is restricted to a few papers and mostly focused onto 
the optical response in both in-plan o 17 i 19 and out-of- 
p^ an9 2 9 4 34 i 3 i 5 i 39 J 40 configurations. 

Electromagnetic eigenmodes in photonic crystal (PC) 
slabs with a periodic pattern have notable differences 
as compared to the ideal reference systems (i.e., not 
waveguide-embedded), which are well known from the 
literature for the cases of both lDM^ and 2Di£ period- 



icities. A most important issue is the light-line problem: 
only photonic modes which lie below the light line of the 
cladding material (or materials, if the waveguide is asym- 
metric) are truly guided and stationary, while those lying 
above the light line in the first Brillouin zone are coupled 
to leaky waveguide modes and are subject to intrinsic 
radiative losses. These quasi-guided modes are actually 
resonances in a region of continuous energy spectrum, 
and for this reason they are more difficult to calculate 
than truly guided modes below the light line. Indeed, 
while the dispersion of guided modes can be obtained 
by a plane-wave expansion with a supercell in the ver- 
tical direction^, the frequencies and especially the losses 
of quasi-guided modes are most commonly studied by 
FDTD calculations Another important feature of pho- 
tonic crystal slabs is the blue shift of the eigenmodes due 
to vertical confinement in the planar waveguide. This 
effect, which is more pronounced for slabs with strong 
out-of-plane refractive index contrast, implies that the 
frequencies of photonic bands and gaps depend on all pa- 
rameters of the planar waveguide (layer thicknesses and 
refractive indices) and can differ substantially from those 
of the reference ID or 2D system. Finally, the eigenmodes 
of photonic crystal slabs can be put in one-to-one corre- 
spondence with those of the reference system only when 
the frequency is sufficiently low for the waveguide to be 
monomode. The cut-off frequency of second- and higher- 
order modes also depends on slab parameters and on the 
photonic lattice. 

In this work we present a systematic study of pho- 
tonic bands, gap maps and diffraction losses for ID pho- 
tonic crystal slabs, that is ID photonic lattices like those 
of a distributed Bragg reflector (see Fig. QJi for the ID 
reference system). These are defined in two types of 
waveguides with strong refractive index contrast: the 
self-standing membrane or air bridge (Fig. ^>) and the 
asymmetric photonic crystal slab in which only the core 
layer is patterned (Fig.[TJ;). The latter structure is typi- 
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(a) Multilayer 




(b) Dielectric membrane 




(c) Asymmetric waveguide 




FIG. 1: Photonic structures studied in this work: (a) Ideal 
one- dimensional photonic crystal, with period a and air frac- 
tion / a i r = h/a. (b) Photonic crystal slab consisting of a 
self-standing, patterned dielectric core (air bridge or mem- 
brane) of thickness d surrounded by air. (c) Photonic crystal 
slab, with the pattern defined in a high-index dielectric core 
of thickness d sandwiched between air and an insulating ox- 
ide substrate. Throughout this work we assume: edict = 12, 

Coxide —2.1, € a ir = 1- 



cally realized with the Silicon-on-Insulator (SOI) system 
but may also be realized with GaAs on an oxide layer. We 
assume the following values of the dielectric constants: 
£dici = 12 for the high-index core layer (as appropriate to 
Si or GaAs below the band gap), e a j r = 1, and e OX ido = 2.1 
(as appropriate for SiC>2 or other oxides). The periodic 
patterning is taken along the x direction and throughout 
this paper we assume k y — 0. The gap maps are calcu- 
lated as a function of air fraction in the core layer and for 
different values of the core thickness, thereby exploring a 
wide range of parameters of experimental interest. The 
dependence of the radiation losses on frequency, polar- 
ization, and air fraction is also calculated and discussed. 

The photonic bands and gap maps of a distributed 
Bragg reflector are obviously well known and are exempli- 
fied in Fig- EI 43 Notice that the photonic bands of Fig.|2t 
(which refer to an air fraction / a ; r = 0.3) as well as the 
gap map of Fig. [!J> are degenerate for transverse elec- 
tric (TE) and transverse magnetic (TM) polarizations 
with respect to the plane of incidence: this degeneracy is 
lifted in a waveguide because the confinement-induced 
shift is polarization-dependent, as was already shown 
experimentally^ One of the goals of the present paper is 
to establish whether a complete band gap for both polar- 
izations can occur in a waveguide-embedded ID photonic 
structure. 



Related concepts have been studied in the context of 
dielectric waveguide gratings, also called resonant grat- 
ing filters. 41 ' 44-59 These kinds of diffraction gratings may 
support guided and leaky modes. The latter are reso- 
nantly coupled to an external light beam and give rise 
to narrow resonances in reflection or transmission, which 
may be used for polarization-dependent filters 4 ^ or for en- 
hanced nonlinear optical effects 4L£& Most of the research 
concentrated on systems with a weak dielectric modu- 
lation, e.g., surface relief gratings for filtering and dis- 
tributed feedback) 41 ^ waveguides with a weak refractive 
index contrast within the core region 45 ^ 4 ^* 5 ^ and/or which 
are modeled by a single Fourier component of the dielec- 
tric function. 44,51-53 ' 58 For an extensive list of previous 
literature along these lines and of the different kinds of 
theoretical methods used we refer to the book by Loewen 
and Popov i 4 ^ In all these cases, which can be treated at 
least qualitatively by coupled-mode theory, the disper- 
sion of the waveguide mode is only weakly modified by 
the dielectric modulation and photonic bandgap effects 
are very small. Specific waveguide grating structures 
with strong refractive index modulation in the plane lead- 
ing to an appreciable photonic gap have been studied 
in Refs. |5j|55j for the case of TE polarization, and in 
Refs. 15711591 for both TE and TM polarizations. In these 
strongly modulated cases a rigorous coupled-wave anal- 
ysis (also called the Fourier modal method) is necessary 
and has been used. We point out that the focus of the 
present work is quite different from all these papers, in 
particular for what concerns the systematic calculation 
of gap maps and losses as a function of frequency and of 
various structure parameters. 

This work is organized as follows. In Sec. II we give a 
short description of the theoretical method for calculat- 
ing photonic bands and intrinsic losses in a waveguide. In 
Sec. Ill we discuss a few examples of photonic mode dis- 




0.5 1 0.2 0.4 0.6 0.8 1.0 

k x a/7t air fraction f air 

FIG. 2: Ideal multilayer: (a) Photonic bands for / a i r = 0.3; 
TE and TM modes are exactly degenerate, (b) Gap map, 
i.e., band edges as a function of air fraction; the value of / a i r 
corresponding to the calculation given in (a) is indicated by 
an arrow. 
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persion, both for symmetric and asymmetric ID photonic 
crystal slabs. Section IV contains detailed results for ID 
gap maps in a membrane or in an asymmetric waveguide 
for different values of the core thickness. In Sec. V we 
present results for intrinsic losses of quasi-guided modes. 
In Sec. VI we give a few closing remarks. 

II. METHOD 

The approac h adop ted here, which was already intro- 
duced in Refs. 1361371 relies on a finite-basis expansion 
in order to transform the second-order equation for the 
magnetic field in a linear eigenvalue problem. The ba- 
sis consists of the guided modes of an effective homoge- 
neous waveguide, where the dielectric constant of each 
layer is defined by the spatial average of the dielectric 
constant e(r) over the photonic pattern. The in-plane 
Bloch vector k is obviously conserved modulo a recipro- 
cal lattice vector G. The off-diagonal components cg,G' 
of the dielectric tensor give rise to band splittings and 
to a folding of the photonic modes in the first Brillouin 
zone. Some of the photonic modes may lie below the 
cladding light line and be truly guided, however most 
(sometimes all) of them fall above the light line in the 
first Brillouin zone. Coupling of these modes to leaky 
modes of the effective waveguide is taken into account by 
time-dependent perturbation theory, which leads to an 
expression for the imaginary part of the mode frequency 
in terms of the photonic density of states at fixed in-plane 
wavevectoriS&SLSI This procedure is formally analogous 
to Fermi's golden rule in quantum mechanics. When 
the waveguide is asymmetric (like in the case of the SOI 
structure), care must be taken to express the leaky modes 
in terms of outgoing states and to relate them to the re- 
spective state densities*^ 

The approximations made in the present approa ch are 
as follows (for a fuller discussion see Refs. I36l37|) . The 
effective dielectric constant of the homogeneous waveg- 
uide, which defines the basis of guided modes for the 
expansion, is chosen to be 

e cff = /air^air + (1 — /airjCdiel , (1) 

for both TE and TM polarizations. This choice is by 
no means unique, although it is the exact definition of 
the effective dielectric constant for TE polarization and 
anyway when the electric field is perpendicular to the di- 
rection of periodicity^ For TM-polarized modes, which 
have electric field components along x and z, the situa- 
tion is more complex^ It is known from the literature 
that the x component of the electric field is subject to an 
effective dielectric constant that is obtained from the in- 
verse averaging rule^ therefore a different choice of e e g 
in the patterned region could be suggested for TM modes. 
Choices of e e g differing from Eq. do not lead to any 
appreciable change of the results above the mode cutoff, 
as we have verified. The frequency position of the cut- 
off does depend on the choice of e e s, especially for large 



air fractions, however a comparison with exact scattering 
matrix calculations 35,36 shows that the average dielectric 
constant defined by Eq. Q gives very good agreement 
with the frequencies and cutoff positions obtained from 
the exact calculations. It should also be noted that the 
electromagnetic field close to mode cutoff is mostly ex- 
tended in the claddings, where the dielectric constants 
are homogeneous for the airbridge and SOI structures 
studied in the present paper. 

The number of reciprocal lattice vectors G is limited by 
a finite cut-off, like for usual plane-wave calculations, 64 
and in addition a restricted number of guided modes of 
the effective waveguide is kept in the expansion. For the 
calculations shown in this work, a number of 31 plane 
waves is usually taken in the basis set and is sufficient 
for convergence with better than percent accuracy. The 
number of guided modes in the expansion is not found 
to be critical in the energy range considered and is usu- 
ally taken to be < 8. For the quasi-guided modes, ne- 
glect of the second-order shift due to coupling to leaky 
modes introduces an error of less than a few percent in 
the photonic frequencies. All these approximations are 
justified a posteriori by the close agreement of the cal- 
culated photonic frequencies with those obtained from 
reflectivity calculations 3 ^ made with the exact scattering- 
matrix method^ Finally, calculating the radiative losses 
of quasi-guided modes by first-order perturbation theory 
is justified by the fact that the imaginary part of the fre- 
quency is much smaller than the real part, as shown by 
the results below. 



III. PHOTONIC BANDS 

The photonic bands of the strong confinement sym- 
metric slab, corresponding to the system schematically 
shown in Fig. ^p, are displayed in Fig. and for a 
core thickness d = 0.4a and an air fraction / a ; r = 0.3. 
The bands are plotted by using dimensionless frequency 
Loa/{2iTc) — a I A as a function of in-plane wave vector 
k x a/TT in the first Brillouin zone. The photonic dispersion 
curves of the patterned structure are compared to those 
of a uniform dielectric slab suspended in air (Figs.|3^ and 
|3Jl) with a spatially averaged dielectric constant given by 
Eq. QJ, that is e e ff = 8.7 in the present case. The 
guided modes of Figs. [5^ and 01 represent the basis set 
for the expansion method discussed in the previous Sec- 
tion. The dispersion for the average dielectric slab is 
presented in the reduced zone scheme, allowing for a di- 
rect comparison with the corresponding photonic bands 
of the patterned waveguide. We have classified the guided 
modes according to mirror symmetry with respect to the 
plane of incidence kz = xz (we use o~ xz to denote this 
operation) and with respect to the xy plane [cr xy opera- 
tion). The modes whose electric field component lies in 
the xy plane are referred to as TE, and are odd with re- 
spect to specular reflection through the plane of incidence 
(o~xz = — 1); the modes whose magnetic field lies in the 
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FIG. 3: Photonic bands for the membrane structure of Fig.^i. The dashed lines represent the dispersions of light in air and in 
the average core layer, (a) TE, (d) TM dispersion curves, folded in the first Brillouin zone, for a uniform dielectric membrane 
with e cf f = 8.7 and thickness d/a = 0.4; (b) TE, (c) TM photonic bands for the patterned structure with / a i r = 0.3, d/a = 0.4. 



xy plane are labelled as TM and are even with respect 
to mirror plane xz (o xz = + 1). 2 « These modes can be 
classified further as even {<j xy = +1) or odd (a xy = — 1) 
with respect to specular reflection through the xy plane, 
thus giving four different types of guided eigenfunctions 
for the electromagnetic field. We can separately compare 
Fig. 0i to|3b and Fig. [3; to|3|i. It is clearly seen that for 
both TE and TM modes the periodic patterning of the 
dielectric slab introduces band gaps around the degener- 
ate points of the average slab dispersion curves (k x — 
and k x = ±7r/a), owing to the off-diagonal components 
of the inverse dielectric tensor. There is one-to-one corre- 
spondence between average slab and ID PC slab modes. 
The first-order modes (TE even and TM odd) have no 
cut-off frequency, as is well known for a symmetric waveg- 
uide. The second-order guided modes have a finite cut- 
off frequency, which is degenerate for TE and TM modes. 
The second-order mode, represented by dotted lines, has 



-1 for TE polarization, while it has a x 



T for 



TM polarization. 



A second point should be remarked by comparing the 
photonic bands of Figs. |3Jd and |2t to the bands of an 
ideal multilayer. The first photonic band gap appears 
between 0.15 and 0.20 in the ideal ID case, and between 
0.22 and 0.28 for the lowest TE mode in the PC slab, due 
to the confinement effect along the vertical (z) direction. 
The gap between the first and the second band opens be- 
tween 0.37 and 0.45 when considering TM modes: these 
values are strongly blue shifted with respect to both the 
multilayer and the TE modes of the PC slab. Thus the 
confinement effect manifests itself in the blue shift of the 
eigenfrequencies of the electromagnetic field with respect 
to the ideal multilayer, and moreover in the removal of 
degeneracy between TE and TM modes: the latter ef- 
fect is due to the stronger confinement of TM compared 
to TE modes in the planar waveguide 4i We also notice 



that all the band gaps, except for the first one, lie in the 
region of guided resonances, and could be experimen- 
tally tested by external reflectance measurements. The 
first band gap, either TE even or TM odd, is instead in 
the region of truly guided modes. A complete band gap 
common to both polarizations can also be seen around 
iva/(2nc) ~ 0.4, where the second-order TE gap overlaps 
the first-order TM gap. As we will see in the next section, 
this is rather a coincidence for ID PC slabs. 

In Fig. 01 we display the photonic bands for the asym- 
metric structure represented in Fig.^J;. The dashed lines 
are the dispersions of photons in air, substrate and ef- 
fective dielectric core. The parameters used in these cal- 
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FIG. 4: Photonic bands for the patterned dielectric-on- 
insulator structure of Fig. ^ with / a ir = 0.3, d/a = 0.4. The 
dashed lines are the light dispersions in the effective core and 
in the upper and lower claddings; a labels the order of the 
guided mode, (a) TE, (b) TM modes. 
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dilations are d/a = 0.4, / a ; r = 0.3, allowing for a direct 
comparison with the results of Figs. |3p and 0;. Owing 
to the asymmetry of the vertical waveguide, a xy is no 
more a symmetry operation: the modes can only be clas- 
sified as odd (TE, Fig. gji) or even (TM, Fig. Hp) with 
respect to the plane of incidence. However, we have indi- 
cated the approximate order of the waveguide mode by 
the index a in Fig. 0] (this can be defined only when the 
modes are well separated in frequency, otherwise mixing 
and anticrossing effects occur). For an asymmetric slab 
there are no modes starting at u> = 04i By comparing 
Figs. Eland 01 we notice that the lowest TE mode of the 
asymmetric ID PC slab is in quantitative agreement with 
the first-order TE mode of the PC membrane; instead, 
the TM modes of the asymmetric slab are somewhat less 
confined than those of the PC membrane. It is important 
to stress that the modes lying between the two claddings 
light lines (oxide and air in this case) are evanescent in 
air, but leaky in the substrate. These modes have finite 
radiative losses, as we will see in section We also 
notice that no complete band gap is present in the asym- 
metric ID PC slab, at variance with the corresponding 
symmetric structure. The results shown in Fig. 01 are 
conceptually similar to Brillouin diagrams calculated for 
TE polarization in the case of an asymmetric waveguide 
grating structured We point out that the described fea- 
tures of photonic band structures for an asymmetric ID 
PC slab were experimentally verified by variable angle 
surface reflectance performed on SOI structures, for both 
TE and TM modest 



IV. GAP MAPS 

In this section we present a complete set of gap maps 
for waveguide-embedded ID photonic crystals. We con- 
sider a band gap as a frequency region where no photonic 
modes exist, either truly guided or quasi-guided above 
the light line. We present the gap maps for modes with 
TE or TM polarizations, i.e., with definite parity with 
respect to the vertical mirror symmetry a xz : this con- 
vention applies to symmetric as well as asymmetric ver- 
tical waveguide structures. For the case of the asymmet- 
ric structure, for which the lowest-order waveguide mode 
has a finite cut-off, only the frequency region above the 
lowest-order cut-off is physically relevant. 



A. Dielectric membrane 

In Fig. [5] we display the calculated gap maps for the air 
bridge structure of Fig. We show the maps for three 
different slab thicknesses, namely d/a — 0.2, 0.4, and 

0. 8. We display in black the true complete band gap, 

1. e., the frequency region in which no photonic modes (or 
resonances) are allowed for any polarization. The gap 
maps are shown for air fraction varying from to 0.7, 
which represent a realistic range for practical realization. 



The solid lines in Fig. [3] represent the cut-off frequency 
of the second-order waveguide mode, which is given by 



uj c a 



2irc 2dy/e c s - e a 



(2) 



and is the same for both polarizations. 

An important feature that we can see from Fig. [5] is 
that for d/a = 0.2, the TE gap map is qualitatively sim- 
ilar to the ideal multilayer one (see Fig. |2fc>) with a blue 
shift arising from the confinement effect. The band gaps 
for TM modes are shifted to much higher frequencies 
and some complete band gaps start to appear only at 
a/A ~ 0.7. The gap map is more complex for d/a — 0.4, 
due to the appearence of higher-order waveguide modes 
at low frequency. Nevertheless, for a/ A < 0.4 the slab 
is still monomode, and a large complete band gap opens 
in a wide range of air fractions (Fig. 0>) . This complete 
gap comes from the overlap of the first TM gap (at the 
Brillouin zone edge, see Fig. EJJ;) and the second TE one 
(at the zone center, Fig. |3p). No complete band gap 
has been found for other values of d/a (calculations not 
shown). For d/a — 0.8 the photonic band structure is 
quite complex because the slab becomes multimode al- 
ready at low frequencies. The band gap in TM modes is 
still present around a/A ~ 0.4, but no complete band gap 
exists because of the presence of second-order TE modes. 
The conclusions from these results are the following: (i) 
the TE gap map in a waveguide resembles the ideal ID 
one only below the cut-off of second-order modes, (ii) 
the TM gap map is very sensitive to the structure pa- 
rameters, and (iii) a complete gap for both polarizations 
is calculated to occur only for a core thickness around 
d/a = 0.4. 



B. Asymmetric waveguide 

In Fig. we show the calculated gap maps of the 
asymmetric ID PC slab of Fig. for core thicknesses 
d/a = 0.2, 0.4, and 0.8. The gap maps of Fig. [5] show 
notable differences as compared to those of the PC mem- 
brane. 

One of the peculiarities of the asymmetric structure is 
the existence of a finite cut-off frequency for the lowest- 
order TE and TM modes. The cut-off frequency as a 
function of air fraction is plotted with dashed lines for TE 
modes, and with solid lines for TM modes. The values for 
the cut-off frequencies obtained by the present approach 
coincide with those following from the expression for a 
uniform, asymmetric slab. The formula 
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(3) 

where r = 1 for TE modes, r = e e ff/e a ir for TM modes, 
and m > is an integer. 

For d/a = 0.2 the asymmetric ID PC slab has only 
first-order TE and TM modes in the whole frequency 
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FIG. 5: Gap maps for the membrane structure of Fig.^j, as 
a function of the air fraction / a ir = h/a. Solid lines represent 
the cut-off frequency of the second-order waveguide mode, (a) 
Core thickness d/a = 0.2, (b) d/a = 0.4, (c) d/a = 0.8. 
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FIG. 6: Gap maps for the asymmetric PC slab structure of 
Fig. 0i. Dashed (solid) lines represent the cut-off frequencies 
of the first- and second-order waveguide modes for TE (TM) 
polarization, (a) Core thickness d/a — 0.2, (b) d/a = 0.4, (c) 
d/a = 0.8. 



range shown. The TE band gaps are again qualitatively 
similar to those of the ideal ID multilayer, with a con- 
finement effect which is close to that of the membrane 
case (Fig. |5^) ; the TM gaps are instead shifted to much 
higher frequencies as compared to the ID multilayer. The 
first TE band gap is in a region below the cut-off of the 
first-order TM mode, thus it may be considered as a com- 
plete band gap. For d/a = 0.4 and 0.8 a second-order TE 
cut-off appears at frequencies around 0.44 and 0.22, re- 
spectively: the TE gap map is similar to that of the ID 
multilayer only below the second-order cut-off frequency. 
The TM gaps are always very different from those of the 
ideal ID case and also quite different from those of the PC 
membrane: TM modes are seen to be extremely sensitive 
to the structure parameters (core thickness and claddings 
dielectric constants). As it can be seen by comparing 
Figs. EJ) andHjjD, the complete band gap for d/a = 0.4 
occurs for the particular case of a ID PC membrane but 
not in the asymmetric PC slab. For d/a = 0.8 a complete 
band gap resulting from the overlap of the first TE and 
TM gaps appears around a/A ~ 0.3 for / a i r > 0.5. 

As a general remark, the numerical results for pho- 
tonic bands and gaps previously shown relate only to the 
real part of the frequency and do not consider the ef- 



fect of coupling to radiative waveguide modes. Thus the 
physical relevance of a photonic band dispersion is ex- 
pected to decrease on increasing the frequency far from 
the light line. We also notice that the concept of mode 
cut-off for resonant modes is not clearly defined when 
radiative broadening is taken into account. For these 
reasons, the gap maps calculated here are expected to be 
more useful in the low frequency region, in particular for 
the band gaps which open below the second-order cut-off 
lines shown in Figs. |S] and [BJ 



V. DIFFRACTION LOSSES 

To complete our analysis of ID PC slabs we have to 
address also the imaginary part of frequency, which gives 
information about the radiative losses due to out-of-plane 
diffraction. This is done by using time-dependent pertur- 
bation theory for the electromagnetic problem, as previ- 
ously discussed in Section [H] We display in Fig. the 
results for parameters d/a = 0.2 and / a j r = 0.3, for both 
a PC membrane and an asymmetric PC slab. In Figs.Et 
and 03 the band diagram and the corresponding imag- 
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FIG. 7: (a) Photonic bands and (b) imaginary part of fre- 
quency for a symmetric ID PC membrane; (c) photonic bands 
and (d) imaginary part of frequency for an asymmetric ID PC 
slab. Parameters are d/a = 0.2, / a i r = 0.3. Solid (dashed) 
lines are TE (TM) modes. The three markers in (b) corre- 
spond to the points marked in Fig. |H| (see text). The arrow 
in (d) denotes the cusp, which corresponds to the second TE 
band in (c) crossing the air light line. 



inary part of frequencies are shown for the symmetric 
ID PC slab. In Fig. [2K the waveguide is monomode for 
both TE and TM polarizations. In Fig. [7jp we show the 
dimensionless quantity Im(w)a/(27rc), corresponding to 
each photonic band of Fig. as a function of mode fre- 
quency. The imaginary part is generally much smaller 
than the real part, indicating the validity of the pcrtur- 
bative treatment adopted. The losses go to zero when 
the mode crosses the light line in air and becomes truly 
guided. It is clear from the figure that the radiative losses 
generally increase on increasing the photonic band index, 
however the behavior of the losses within a given photonic 
band is nontrivial and has to be studied in each specific 
case. The guided resonances at the Brillouin zone center 
present a quite interesting behavior. In fact, the second 
and the fourth TE bands have zero linewidth at k x = 
(their frequencies are uja/(2nc) — 0.45 and 0.81, respec- 
tively) while the third and the fifth band have finite ra- 
diative widths at k x = (wa/(27rc) = 0.56 and 0.94). 



A similar behavior holds also for TM modes. These nu- 
merical results could be probed by variable angle surface 
reflectance experiments made on ID PC membranes: the 
imaginary part of the frequency can be extracted from 
the linewidth of spectral structures in reflectance that 
correspond to photonic modes»i^2fi 

In Figs.0: and0I we show the results for an asymmet- 
ric ID PC slab with the same thickness and air fraction. 65 
The radiative losses shown in Fig. Qi display quite the 
same behavior as in the membrane case. They are about 
two times larger than the corresponding losses of Fig. [3d: 
this is due to the asymmetry of the vertical waveguide, 
which implies that a quasi-guided mode above the light 
line is coupled to all radiative modes of the effective 
waveguide at the same frequency, without the parity se- 
lection rule which holds instead in the symmetric case. A 
similar behavior was found in the calculation of spectral 
properties of deep ID gratings^ Moreover, the higher- 
order modes (either TE or TM) now have a finite Im(w) 
also at k x = 0: this is due to the additional diffrac- 
tion channels for radiative losses which are present in 
the asymmetric waveguide. Moreover, we notice that the 
modes whose frequencies lie between the light lines of air 
and oxide claddings are not truly guided, i.e., they are 
evanescent in air but leaky in the substrate. Thus, the 
crossing between a band and the light line in air does 
not cause the linewidth of the photonic resonance to go 
to zero: rather, Im(o;) has a cusp (marked by an arrow in 
Fig. 01) when the light line in air is crossed. Similar fea- 
tures can by recognized in Fig. up at higher frequencies: 
they arise whenever a photonic mode crosses a cladding 
light line folded in the first Brillouin zone. These notable 
features of Im(w) are not a numerical artifact, but rather 
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FIG. 8: Imaginary part of photonic frequencies as a func- 
tion of the air fraction for the PC membrane of thickness 
d/a — 0.2; the three curves correspond to different wavevec- 
tors (k x — 0, k x = n/2a, k x = ir/a) in the first Brillouin 
zone of the third TE photonic band. The three points marked 
on the curves correspond to the three markers of Fig. |7|d, 
where fair = 0.3. 
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they correspond physically to the opening or closing of 
diffraction channels for radiative losses. 

We also found that the imaginary part of frequency in- 
creases on increasing the air fraction in the investigated 
range, as shown in Fig. [SJ where Im(o;a/27rc) is plot- 
ted for the PC membrane of thickness d/a — 0.2. A 
similar behavior is found also for the asymmetric struc- 
ture (not shown here) and an increase of the photonic 
mode linewidth with the air fraction was already stated 
experimentally. 21 - 35 The three curves of Fig.[S]correspond 
to the evolution of the losses for the third TE photonic 
band at three different points in the first Brillouin zone. 
When f a i r — 0.3 the corresponding photonic band has 
frequencies around wa/(27rc) = 0.6 (see Fig.EJi), an d the 
three points marked in Fig.|SJcorrespond to those marked 
in Fig. [TJd. Notice that the losses vary in a logarithmic 
scale and become extremely small either towards the ho- 
mogeneous waveguide limit at low air fraction or close to 
the Brillouin zone edge. While the filling fraction depen- 
dence of the losses is similar for all bands and polariza- 
tions, the wavevector dependence changes from band to 
band, as it appears from Figs.[7t>,d. It can be concluded 
that for the present waveguide-embedded ID photonic 
structures, the spectral linewidth of quasi-guided pho- 
tonic modes can vary by several orders of magnitude and 
it depends in a nontrivial way on the structure parame- 
ters as well as on the angle of incidence, mode index, and 
polarization. 

VI. CONCLUSIONS 

We have studied the photonic bands and the gap maps 
of one-dimensional photonic crystal slabs made of a high 
refractive index material sandwiched between low index 
claddings, thus providing strong confinement of electro- 
magnetic waves along the vertical direction. The method 
adopted yields the frequencies of photonic modes both 
below and above the light line, and therefore allows treat- 
ing guided and quasi-guided modes on the same footing. 

Photonic modes in a ID photonic crystal waveguide 
can be put in one-to-one correspondence with those of 



the ideal ID reference system only when the waveguide 
is monomode. Under this condition, the TE gap maps 
are qualitatively similar to the ideal ID ones consider- 
ing the blue shift in the waveguide. This confinement 
effect is considerably more pronounced for TM-polarizcd 
than for TE-polarized modes, thereby leading to a po- 
larization splitting which depends in a sensitive way on 
the structure parameters. As a consequence, a complete 
band gap common to TE and TM polarizations is gen- 
erally not found in ID photonic crystal slabs, except for 
special values of the parameters. 

The radiative losses of guided resonances due to out-of- 
plane diffraction depend in a sensitive way on waveguide 
parameters, mode index, and frequency. For some modes 
the imaginary part of the frequency vanishes at the Bril- 
louin zone center, thus even above the light line it is 
possible to find photonic modes with very low losses. In 
general, the losses are predicted to be higher for Silicon- 
on-Insulator structures as compared to photonic crystal 
membranes, due to the asymmetry of the planar waveg- 
uide. All these results are related to the presence of 
diffraction channels for radiative couplings. 

The quantitative results for the complex frequency dis- 
persion of quasi-guided modes may be experimentally 
tested by performing reflectance or transmittance mea- 
surements with light incident on the PC slab surface. 
Modes below the light line can be probed in waveguide 
transmission experiments. Moreover, the results of the 
present work may be useful for designing ID PC slabs 
with gaps at specified frequencies, regions of monomode 
propagation, and low diffraction losses. They might also 
be used for the design of advanced ID structures like fil- 
ters and microcavities. The same theoretical approach 
can be applied to PC slabs with various 2D patterns. 
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